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Abstract We study Anderson localization of single particles in contin- 
uous, correlated, one-dimensional disordered potentials. We show that 
tailored correlations can completely change the energy-dependence of 
the localization length. By considering two suitable models of disorder, 
we explicitly show that disorder correlations can lead to a nonmono- 
tonic behavior of the localization length versus energy. Numerical cal- 
culations performed within the transfer-matrix approach and analytical 
calculations performed within the phase formalism up to order three 
show excellent agreement and demonstrate the effect. We finally show 
how the nonmonotonic behavior of the localization length with energy 
can be observed using expanding ultracold-atom gases. 

1 Introduction 

Coherent transport in disordered media shows considerable interest in condensed- 
matter physics, with applications to normal solids [1,2], superconductors [3] and su- 
perfluids [4]. Coherent process can lead to the spatial localization of wave functions 
as a result of a subtle interference effect between multiple scattering paths, which 
survives disorder averaging. This effect, known as Anderson localization (AL), was 
first predicted for electronic matter waves [5]. It was later shown to be a universal 
phenomenon in wave physics [6] , which permitted the first evidence of AL of classi- 
cal waves [7,8,9,10,11]. The observation of AL in ultracold gases in one (ID) [12,13] 
and three (3D) [14,15] dimensions has triggered a renewed interest on matter wave 
localization, and paves the way to further investigation of many open questions [16]. 

Correlated disorder makes AL fascinatingly rich. Disorder correlations can change 
the localization properties, not only quantitatively but also qualitatively. For instance, 
disorder correlations with a finite support in momentum space were shown to induce 
effective mobility edges in ID disorder [17,18,19], which was used to create materials 
with alternating localizing and almost transparent frequency windows [20] , to enhance 
localization [21] in microwave systems, and to propose realization of atomic band-pass 
filters [22] . Such correlations are also responsible for algebraic localization of matter 
waves with broad energy distributions [12,18,23]. In addition, tailoring the disorder 
correlations can lead to a counter-intuitive behavior of the energy dependence of the 
localization strength [24] . It can serve to discriminate quantum versus classical local- 
ization, which is of particular interest for ultracold atoms where disordered potentials 
can be controlled [25,26]. 
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The latter effect is the subject of the present paper. For the sake of clarity, con- 
sider matter waves in free space and subjected to a ID disordered potential. The 
macroscopic behavior of AL is intimately related to the microscopic properties of 
single scattering from the asperities of the disorder. This shows up in the strong 
dependence of the localization length L loc (E) - i.e. the characteristic length scale 
of exponential decay of localized wavefunctions - on the energy E and the Fourier 
component of the disordered potential V at twice the typical particle wave vector, 
k E = yj2mE/h: The ID localization of a particle is dominated by the interference of 
quantum paths that are backscattered twice in the disorder with the forward prop- 
agating one. For weak disorder, the leading term is L loc (E)^ 1 oc C2(2k E )/E, where 
C2 is the structure factor (correlation function of the disordered potential in Fourier 
space, see below) [27]. For most models of disorder (e.g. for (^-correlated, Gaussian- 
correlated or usual speckle [25] disorder), C2 is a constant or decreasing function, so 
that L loc (E) is an increasing function of the particle energy E, which finds an intuitive 
interpretation as the higher energy the weaker localization. 

Structured disorder correlations can however completely change the picture. If the 
disordered potential exhibits structures on a length scale of the order of fc" 1 , the scat- 
tering might not strictly increase with energy 1 , and L loc {E) can counter- intuitively 
decrease for increasing E [24]. In this work, we study this effect for particles in con- 
tinuous disordered potentials with tailored correlations in one dimension. We first 
calculate the localization of the single-particle eigenstates, using both numerical cal- 
culations, based on the transfer-matrix approach, and the so-called phase formalism, 
which is well suited for perturbative expansion in ID transmission schemes for a 
matterwave of fixed energy. Lowest-order analytical calculations reproduce the main 
physics and allow design of disorder correlations to realize the desired effect [24] . We 
study two alternative possibilities, which require only slight modifications of exist- 
ing experimental schemes [12,29,30,31]. These calculations however show significant 
deviations with numerical data, but we show that they are quantitatively accounted 
for by next-order calculations. We finally discuss how to observe the nonmonotonic 
behavior of the localization length with energy with expanding ultracold-atom gases, 
and explicitly show that standard schemes should be adapted. 



2 Tailoring disorder correlations in speckle potentials 

Disorder realized by optical speckle [32,25,26] is relatively easy to obtain and its sta- 
tistical properties are -to some extent- controllable, hence making a good candidate 
for tailoring correlations. Such disordered potentials have been used in several experi- 
ments to investigate the effect of disorder in ultracold-atom systems [12,29,30,31,14,15]. 
A ID speckle pattern [32] is obtained by shining a coherent light beam onto an elon- 
gated aperture with a ground-glass plate diffuser, and observing the diffraction pat- 
tern that is created, for example in the focal plane of a lens (see Fig. 1). The role of 
the diffuser is mainly to imprint a random phase on the electric field at each point 
of its aperture. It follows from direct application of the laws of diffraction, that the 
electric field at a point z on the focal plane, E (z), is the sum of many complex in- 
dependent random variables that corresponds to the components originating from 
every point of the plate and interfering in z. It is then a complex Gaussian random 
variable. In ultracold atom experiments, the atoms are sensitive to the light intensity, 
proportional to |£(z)| 2 . More precisely the atoms experience a potential that can be 



1 This effect is similar to the well-known example of single scattering from double barriers 
(see for instance Ref. [28]) 
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Figure 1. Optical scheme to create a 
speckle pattern: A laser beam is diffracted 
by a ground-glass plate diffuser (D) with 
a ID slit of pupil function I (Z), where 
Z spans the diffuser. The latter imprints 
a random phase on the various light paths. 
The intensity field, I(z) , observed in the fo- 
cal plane of a converging lens, is a speckle 
pattern, which creates a disordered poten- 
tial V(z) for the atoms. 



written V(z) = V R \£(z)\ 2 /\£\ 2 - 1 , so that V = and V 2 = V£ [25,26]. The ampli- 
tude of the disorder, V R , can be positive (so called 'blue detuned speckle') or negative 
('red detuned speckle'), depending on the sign of the difference of the atom and laser 
frequencies. More generally, all statistical properties of the disordered potential V{z) 
follow from basic laws of optics. For instance, one can show that the Fourier transform 
of the autocorrelation function of the potential is the auto-convolution of the pupil 
function (intensity pattern exiting the diffuser): 

C 2 (k) oc fdZlJz- ^k) ljz+ ^k) , (1) 



4tt ) \ 4tt 

where A is the laser wavelength and I is the focal length [32]. For a thin slit of 
length 2R uniformly lit, which is the usual manner to obtain a speckle, we have 
I D (Z) = I O(R - \Z\), with the Heaviside function \0(x) = 1 if x > and 
otherwise], and we find 2 

C2(k)=7rV£a R [l-\k\a R /2] m , (2) 

where [x]^ = xO(x) and <r R = X 1/2t:R is the disorder correlation length. Therefore, 
the power spectrum of the speckle pattern is strictly decreasing with k. Note that this 
model has a cut-off in the Fourier components at k c = 2a~ 1 . This is a characteristics 
of speckle potentials, which in particular leads to the existence of effective mobility 
edges, as discussed in Refs. [18,19,33,34]. 

Modifying the pupil function I (Z) (i.e. changing the aperture of the diffusive plate 
or the spatial profile of the incident beam) allows us to tailor the disorder correla- 
tions [22,24]. In this work, we consider two configurations. In the first configuration, we 
propose to put a mask of width 2r at the center of the aperture, creating a double-slit. 
When doing so, a gap is created in the pupil function [I D (Z) = IoO(\Z\ — r)0(R—\Z\)], 
leading, for p = r/R > 0, to an increase of the integral (1) on a certain interval of k, 
which is all the more marked that p is large: 

C^ k ) = ^-P- l^'^e + \ t 1 - P ~ 11*1^ ~(P+ 1)11© }• ( 3 ) 

In the second configuration, we propose to illuminate an infinitely-long slit by two 
mutually coherent Gaussian laser beams of waist w along Z and centered at Z = 
±A/2. The two-point correlation function is then given by 



(fcj R -K ) 2 (fc<T R ) 2 (fcg R + lt Q ) 2 

e * +2eJ+e 4 



(4) 



We use the conventions /(k) = J f(u)e ZKU du. 
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with cr R = X 1/ttw and n = 2A/w. For k > 3.7, this function also increases on a 
certain interval of k. We will see in the following that the increase of G^k) in both 
the double-slit and the double- Gaussian configurations can lead to an enhancement 
of localization with energy. 



3 Anderson localization in ID tailored speckle potentials 

We now study AL of a single particle in these disordered potentials using the so- 
called phase formalism [27], which allows for efficient perturbative expansions [33]. 
Consider a particle of mass m and given energy E in the disordered potential. The 
corresponding eigenstate 4> E (z) is governed by the ID Schrodinger equation 

- (h 2 /2m)d 2 z Mz) + V(z)<f> E (z) = E<t> B (z). (5) 

Then, write the (real- valued) eigenfunction 4> E (z) = r(z) sin [6(z)] and its spatial 
derivative d z (j) E (z) = k E r(z) cos [0(z)]. In this representation, Eq. (5) transforms into 
the set of equations 

d z 9{z) = k E (l - S m 2 [9(z)]) (6a) 



In 



dz'^ff-sin[2e(z% (6b) 



"(0)j , 

Equation (6a) can be solved in the form of a Born-like perturbative series for the 
phase 9 in powers of the external potential V . Then, the Lyapunov exponent j(E) = 
L loc (E)^ 1 = lim^i^oo ln[r(z)]/|z| can be calculated at each order in V by inserting 
the result of Eq. (6a) into Eq. (6b). It yields 3 j(E) = E«> 2 7 ( ™ ) (^) with 



h{K) = i^3l, (8) 



each function f n depends on the n-point correlation function of the disorder, C n (z\, z n -x) = 
V(Q)V(zi)...V(z n -i) [33]. In particular, the leading term of the series (Born approx- 
imation) is 

I C 2 (2 K /q R ) 

This term generally captures most localization properties in ID disorder, in par- 
ticular the effect of the tailored correlations we consider here. Previous work has 
however shown some discrepancies between analytic calculations in the Born approx- 
imation and numerical calculations [24]. Therefore, we will also include the next-order 
term [33] 

M") = 4f d - f dv C ^ W ^ sin M . (9) 

In addition, we will perform numerical calculations for the transmission through the 
ID disordered potential of a particle governed by Eq. (5), using transfer matrix tech- 
niques [35]. 

Let us start with the standard (single-slit) configuration where the pupil func- 
tion is uniform and nonzero in the interval —R<Z< +R, which corresponds to 
the correlation function given by Eq. (2). Figure 2 shows the Lyapunov exponent in 
this configuration, for both blue- and red-detuned speckle potentials, with parame- 
ters relevant for current experiments [12]. The numerical data (blue squares and red 



Note that the first-order term vanishes because V = 0. 
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Figure 2. (Color online) Lyapunov expo- 
nent for both blue- and red-detuned speck- 
les with the standard correlation func- 
tion (2) and V R = ±0. 01625(/l 2 /mcr|). 
Shown are the numerical results extracted 
from a transfer matrix method (space step 
of Az — 0.1<7r and total system size of 
Ltot = 411775<7 R with random initial con- 
ditions) averaged over 5000 disorder real- 
izations (blue squares: V R > 0, red dia- 
monds: Vr < 0), and analytical results ob- 
tained from the phase formalism up to or- 
der 2 (solid black line) and up to order 3 
(dotted blue and red lines). 

Figure 3. (Color online) Probability dis- 
tributions of the Lyapunov exponents ob- 
tained by transfer matrix calculations (dot- 
ted lines) with the disorder parameters of 
Fig. 2, Vr > (space step of Az — O.Ictr 
and total system size of Ltot = 13333<j r 
with random initial conditions), and vari- 
ous energies (indicated on the figure). The 
solid lines are the corresponding theoretical 
distributions (10). 



diamonds) are averaged over 5000 realizations for each value of k E . The analytic calcu- 
lations of the Lyapunov exponent in the Born approximation [y^ (E) given by Eq. (7) 
with n = 2 and Eq. (8); solid black line], which do not depend on the sign of Vr, fairly 
reproduce the numerical data. In this standard configuration, both numerics and ana- 
lytic calculation in the Born approximation confirm that the Lyapunov exponent, i.e. 
the localization strength, decreases with increasing particle energy, hence following 
the intuitive behavior. We however find a significant discrepancy, which depends on 
the sign of Vr between, the numerics and the analytics. A similar discrepancy was also 
observed in Ref. [24] but the numerics were not calculating 7 (£7) in a direct way as in 
the present work. As is seen on Fig. 2, this discrepancy is very well accounted for by 
analytic calculations to the next order in the perturbative series, ^ 2 \E) + j( 3 \E), 
where 7^ 3 ^(£') is given by Eq. (7) with n = 3 and Eq. (14) with p = 0. Note that 7' 3 ) 
depends on the sign of Vr. 

Figure 2 shows the average values of the Lyapunov exponent, which we write j(E) 
in this paragraph. For the sake of completeness, we show in Fig. 3 some distributions 
of the Lyapunov exponents j(E) obtained for the different realizations of the potential 
at four values of the energy, together with the Gaussian distribution 

F(7) = 7ib; exp 



(7-7)' 
2AI 



(10) 



with Ay = y/j/Ltot an d Ltot the system size, which is expected for ^-correlated 
disorder when L to t I/7. For low energy, we find a good agreement. Since -f(E) = 
— ln(T)/2L t ot, where T is the transmission probability of the wave through a sample 
of finite length Ltot, it shows that the transmission probability T follows a log- normal 
distribution also in the correlated disorder we are considering. At high energy (k E a n > 
0.8) the numerical distributions differ from Eq. (10). This is expected because for 
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Figure 4. (Color online) Lyapunov expo- 
nent for speckles with the tailored correla- 
tion function (3) with p = 1/3, using the 
same methods and parameters as in Fig. 2. 
Blue squares: numerics with Vr > 0; Red 
diamonds: numerics with Vr < 0; Solid 
black line: order 2 in the phase formalism; 
Dotted blue and red lines: order 3 in the 
phase formalism. 

Figure 5. (Color online) Lyapunov expo- 
nent for speckles with the tailored correla- 
tion function (4) with k = 8.88 and Vr = 
±0.72(h 2 /mai), using the same methods 
as in Figs. 2 and 4. Blue squares: numer- 
ics with Vr > 0; Red diamonds: numerics 
with Vr < 0; Solid black line: order 2 in 
the phase formalism; Dotted blue and red 
lines: order 3 in the phase formalism. 



fc E (T H = 0.81, we find I/7 ~ 5000(T R , which is of the order of magnitude of L to t = 
13333(7 R . 

As discussed above, for standard disorder, i.e. with a power spectrum C2 that is a 
constant or decreasing function of k, the Lyapunov exponent ^ 2 ^(E) decreases with 
the energy E. Let us now consider the double-slit configuration, of correlation function 
given by Eq. (3). Inserting the latter into Eq. (8), we find that, for p > 0.25, 7 (2) (#) 
shows an increase in a certain interval of k, which is all the more pronounced that 
p approaches 1. This indicates that in those tailored potentials, j(E) can counter- 
intuitively increase with energy. In order to study this effect precisely, we performed 
numerical (transfer-matrix approach) and analytical (phase formalism approach) cal- 
culations for the considered tailored speckle potential with p = 1/3, as done for the 
standard speckle potential. The results of the numerical calculations (blue squares 
and red diamonds) and of analytical calculations up to order three in the phase for- 
malism (dotted blue and red lines) are shown on Fig. 4. They confirm that, for both 
blue and red detunings, j(E) exhibits an increase with E for fc £ cr R £ [0.35,0.6]. As 
for the standard speckle potential, the numerical results follow the trend of the 2 nd 
order term in the phase formalism (solid black line in Fig. 4). The 3 rd order term 
is given in the appendix [Eq. (14)]. For k E a R € [0.4, 1], it accounts very well for the 
discrepancy between the numerics and the 2 nd order term (solid black line). For low 
energy (i.e. k E <r n < 0.4), the 4 th and higher order terms play a more important role, 
which is expected as we are approaching the limits of validity of the perturbativc 
development. 

In the above double-slit configuration, j(E) has a slope break near its maximum 
[see Eqs. (3), and (14) and Fig. 4], which is reminiscent of the sharp edges of the pupil 
function. As it will presumably be inconvenient for experimental observations, we now 
consider the double-Gaussian configuration, which is obtained using two mutually 
coherent Gaussian beams shone onto an infinite diffusive plate, giving the power 
spectrum (4) [24]. For this configuration, ^ 2 \E) shows an increase when n > 5.3, 
which is all the more marked than k is large. In Fig. 5, we show the Lyapunov 
exponents obtained in this case for n„ = 8.88, with transfer matrices (blue squares 
and red diamonds) and with the phase formalism, up to order 2 (solid black line) and 
up to order 3 (dotted blue and red lines; see Eq. (15) in the appendix for the 3 rd order 
term). In this configuration, we recover the same trend as in the other configuration, 
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namely the Lyapunov exponent shows a significant increase in a given energy window, 
k E a n £ [2.3,4.2], the second order term, j( 2 \E), captures the main physics, and the 
discrepancy between the numerical results and (E) are well accounted by the third 
order term, except at very low energy where the perturbative expansion breaks down. 
As expected, the behavior of 'y(E) is smoother for the double-Gaussian configuration 
compared to the double-slit configuration. 



4 Observation schemes with ultracold atoms 

In order to probe the nonmonotonous behavior of *y(E) discussed above, one can use 
ultracold atoms, which proved a good means to observe ID AL of matter waves with 
pseudo-periodic [13] and speckle [12] potentials. The preceding calculations of the 
Lyapunov exponent (a self-averaging quantity) directly apply to a ID transmission 
scheme of a wave with fixed energy E. In ultracold-atom experiments [12,13,14,15] 
however, a matter-wavepacket with a broad energy distribution should be considered, 
and the measured quantity is the density profile obtained after releasing the atoms 
in the disorder, which is not directly related to the above calculations. The average 
stationary density of a noninteracting atomic gas, with initial negligible width, after 
evolution in the disorder reads [18,19,23] n 0o (z) = J dE r D E (E)P 00 (z\E), where T> E (E) 
is the energy distribution of the atoms and 



7T 2 7 f°° 

P ca (z\E) = / du u sinh(7ru) 

8 Jo 



1 + u A 



1 + cosh(7ra) 

xcxp{-(l + U 2 ) 7 |z|/2}, (11) 

with 7 = j( 2 )(E) given by Eqs. (7) and (8), is the probability of quantum diffusion 
calculated in the weak disorder approximation [36,37,38]. 

A first attempt to observe the nonmonotonous behavior of 7 (i£) may be to consider 
the experimental scheme of Ref. [12]. In this case an interacting condensate is first 
produced in the Thomas- Fermi regime in a harmonic trap of frequency w, and the trap 
is then switched-off at time t = 0. In a first stage the expansion of the atoms is driven 
by their interaction energy, and one can neglect the disordered potential. For t 3> 1/w, 
it produces an almost noninteracting gas with momentum distribution [39,40] 

A(p) = (3/4p cut )[l-bM ut ) 2 ] e - (12) 

We have performed numerical integration of the time-dependent Schrodingcr equation 
for a particle with the initial momentum distribution (12) in the two-Gaussian tailored 
disordered potential with correlation function (4) and disorder parameters as in Fig. 5, 
for six realizations of the disordered potential [three with blue (V R > 0) and three 
with red (V R < 0) detuning] 4 . After averaging the stationary density profiles over 
the six realizations, we fit In [Poo (z)] as given by Eq. (11) to lnfnT^z)] with 7 as the 
only fitting parameter 5 . The results, plotted on Fig. 6, show that the fitted Lyapunov 
exponent (black dots) slightly decreases with k cut and saturates roughly beyond the 
minimum of the calculated 7^. This is because the long distance behavior of n 00 (z) is 
dominated by the energy components with the largest localization lengths, i.e. those 



4 As in Refs. [23,24], we use a Crank-Nicolson algorithm with numerical parameters: space 
step Az = 0.2a R , time step At = 0.025ft/ E c , boxes of size 80 x 10 3 ct r . 

The fits are performed in the space windows — 300<t r < z < — 50<r R and +50<r R < z < 
+300<J R , corresponding to an experimentally accessible width of 1 mm for cr R = 1.6 \x.m. 
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Figure 6. (Color online) Lyapunov ex- 
ponent versus particle energy for speckles 
in the double-Gaussian configuration with 
« = 8.88 and V R = ±0.72(ft 2 /mo£), ob- 
tained from Eq. (8) (solid black line) and 
from fits of Eq. (11) to the numerical data 
(points). The figure shows the values of y/u 
extracted from the average profile rioo(z), 
with initial Thomas-Fermi [TF, Eq. (12), 
full dots] and Gaussian [Eq. (13), open cir- 
cles] momentum distributions. 



with the smallest 7(E) [18,19,23]. Therefore this scheme does not enable us to probe 
the region where "f(E) increases. 

In order to observe the upturn of j(E), one can use an atomic energy distribu- 
tion much narrower in energy and strongly peaked at a given E at , so that 7100(2) ~ 
P oa (z\E a- t)- As discussed in Ref. [24], it can be realized by either giving a momentum 
kick to a noninteracting initially trapped gas or using an atom laser, both with a 
narrow energy width [41,42,43]. Using the momentum distribution 

A(p) = (l/\/2^w) cxp [-(p - p at f/2pl] . (13) 

with parameters relevant to current experiments, one can then extract the values of 
7 (.Eat) by the same fitting procedure as above. The results, displayed on Fig. 6 (grey 
circles, reproduced from Fig. 2(a) of Ref. [44]), show a strong increase of j(E) as a 
function of fc at , which follows quite well Eq. (8). The scheme, which require a small 
change in current experiments hence allows one to directly observe the nonmonotonic 
behavior of 7(E) induced by the tailored correlations. 



5 Conclusion 

In summary, we have studied Anderson localization of noninteracting ultracold mat- 
terwavcs in correlated disordered potentials. Tailoring of the disorder correlations can 
strongly affect the localization behavior, namely the localization length can increase 
with particle energy. Here we focused on the ID case, providing more details compared 
to earlier works [22,24]. In particular, we show that although the Born approxima- 
tion reproduces the overall behavior of the localization length, next-order terms are 
significant. We have found good agreement between exact numerical calculations and 
perturbation theory up to order three. We have compared two disorder configurations 
('double-slit' and 'double-Gaussian') to realize the above effect. We finally discussed 
how to observe the effect with expanding ultracold atoms, explicitly showing that the 
scheme usually used in experiments [12,13,14,15] should be adapted. 
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de Recherche sur les Atomes Froids (IFRAF). We acknowledge the use of the computing fa- 
cility cluster GMPCS of the LUMAT federation (FR LUMAT 2764). 
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Third-order term in phase-formalism calculations For the double-slit configuration, 
the / 3 function that intervenes in Eq. (7) reads 



/s(«) 



8(1 -Pf 



0(t £ -'O/3,i(k) + e(Y-|»- / r 1 l) 







(k - ^hA*) + e(^-K) / 3 , 3 («) 



(14) 



with 



/ 3 ,i(k) - 5 (2) - fl[2(l - «)] - 2. 9 (2k) + 2g (1 - p) - g (2p) -2g(l-p- 2k) 
+g (1 - 2k + p) + .g [2 (« + p)] - g (1 + 2k + p) , 



and 



/s,a(K) = -ff(2/c) - <? (1 - p) + <? (2p) + .9 (2k + 1 - p) 
+.9 (2k - 2p) - g (1 + p) + g (l-2n + p), 



/ 3 , 3 (k) = 5 (2) - g[2(l - k)] - g(2n) +g (1 - p) - g (-1 + 2k - p) 
-g{l + p)+g{2n-l- p) 



where g(x) = xhi(x). 

For the double-Gaussian configuration, we find 



/a(«) 



1/1 



47re fe i "b crfi 



8^ \2fciO-R 

1 (K + 2fc c <T R ) 2 / C\ T 

~3 / K + ZK c tJ R 



+ 27re fc i "r x 



(15) 



t r erfi 

2Kfc, 



1 (n-2k c v- R )^ 

3 fcTZ3 j2 

e i r era 



k — 2£) c <7 

\/3fc, 



cr R / J 



^ftfc ( 1 (K+fc c <T R ) 

+2e fc r cosh( — — -)< e 3 fc i"R crfi 
fcfcr R I 



k + k c a R 
V3fc,cr R 



+e 3 fc i "r erfi 



K k c <7 Ti 



-)} 

CTr / J 



where erfi (2) = —i erf (iz) is the imaginary error function, with erf (x) = 2 J* dt e 1 /t/tt 
the error function. 
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